On negative streamers: a deterministic approach 
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In dielectric breakdown, the phenomena of streamers formation and spontaneous branching is 
commonly observed. A deterministic negative streamer model is reviewed in this article. We recently 
have proposed that this reduced model is able to explain the branching phenomena, due to a 
Laplacian instability quite like in viscous fingering. Particular attention is paid to the physical 
interpretation of the model. 
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I. INTRODUCTION 

Atoms of any substance are electrically neutral. If we 
apply an electric field to a volume filled with neutral par- 
ticles, the electric current will not flow through that vol- 
ume. Since no charged particles are present, there will 
not be any electric current, which is in fact the directed 
motion of charged particles. Thus, the volume, filled with 
atoms of any substance, is an almost ideal insulator. 

The air is a good example of such an insulator. Ev- 
ery cubic centimetre of air contains roughly speaking 
2.7 X 10^^ molecules of oxygen (O2), nitrogen (N2), 
vapour (H2O) and some other gases. Every atom of oxy- 
gen contains 8 positively charged protons and the same 
number of negatively charged electrons. Every atom of 
nitrogen contains 7 protons and 7 electrons. It might 
seem there are more than enough charged particles, but 
those particles are bound by powerful electric forces to 
form electrically neutral atoms and molecules, and as a 
result of which the air is a perfect insulator. 

However, if a strong electric field is applied to matter 
of low conductivity and some electrons or ions are cre- 
ated by some agent, then these few mobile charges can 
generate an avalanche of more charges by impact ioniza- 
tion. A low temperature plasma is being created, result- 
ing in an electric discharge. Examples range from natu- 
ral phenomena like the familiar lightning and St. Elmo's 
fire to lamps (neon tubes, hight brightness flat computer 
and TV screens) and industrial plasma reactors for com- 
bustion gas cleaning, ozone generation, etci There have 
been a huge development of technologies based on this 
phenomena. Laser pumping and ion chambers used for 
particle detectors arc other examples one can find. 

On the other hand, the understanding of the basic 
mechanisms of the electric discharge is a challenging 
problem where ideas from nonequilibrium thermodynam- 
ics, atomic physics, electromagnetism and pattern forma- 
tion come into play. Discharges are nonequilibrium pro- 
cesses occurring in initially nonionized matter exposed 
to strong electric field. Depending on the spatiotcmpo- 
ral characteristics of the electric field and on the ion- 
ization and charge transport properties of the medium, 
discharges can assume many different modes of appear- 
ance. Phenomenologically discharges can be classified in 



stationary ones, such as arc, glow or dark discharges and 
transient ones, such as leaders, initial stages of sparks 
and streamers.^ As a warning, the distinction between 
the various discharge phenomena seems to vary among 
authors. 

A streamer is a sharp nonlinear ionization wave prop- 
agating into a nonionized gas, leaving a nonequilibrium 
plasma behind. They have been reported to appear 
in early stages of atmospheric discharges^ like sparks 
or sprite discharges 1^ It is observed that streamers can 
branch spontaneously, but how this branching is pre- 
cisely determined by the underlying physics is essen- 
tially an open question. We have proposed in recent 
work^ a branching mechanism which is qualitatively dif- 
ferent from other ideas. The older concept of dielec- 
tric breakdown can be traced back to Raether's worki& 
In his model he introduced the rare long ranged photo- 
ionization events. Some stochastic models for dielectric 
breakdown have been proposed and studied since theni 
We used a fully deterministic model with pure impact 
ionization and it becomes a surprise streamers get un- 
stable and develop branching. The mechanism for this 
branching is related to a Laplacian interfacial instability. 

In this paper we will start introducing a basic fluid 
model incorporating the physical ingredients to describe 
a nonattaching gas such as nitrogen under normal con- 
ditions. In Sec. Ill some numerical simulations based on 
this model are shown. In Sec. IV and the ones to fol- 
low, the streamers branching is investigated by analyti- 
cal means. We summarise the main results for stationary 
planar fronts and in Sec. V for shock fronts. In Sec. VI 
we set up the framework of the linear perturbation anal- 
ysis for transversal Fourier modes, first the equation of 
motion and then the boundary conditions and the so- 
lution strategy. In Sec. VII we obtain the asymptotic 
behaviour of the dispersion relation and finally, we end 
with a summary and prospect of future work. 



II. THE STREAMER MODEL 

In this section we present a minimal streamer model, 
i.e., a "fluid approximation" with local field-dependent 
impact ionization reaction in a nonattaching gas like ar- 
gon or nitrogeni^ It treats the dynamics of the free elec- 
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trons and positive ions in a homogeneous gas at rest. In 
detail, this is as follows: 

(i) The ionization by electrons in the gas is essentially 
the primary process in all spark discharge. Initially, an 
electron liberated by any outside agents, as e.g radiation, 
is accelerated in a strong local field. It collides into a neu- 
tral molecule and ionises it. The result is a generation 
of new free electrons and a positive ion. In general, this 
process is determined by the rate of gain of energy of 
the electrons and their ability to produce ionization once 
they have sufficient energy. The energy gain depends on 
the field strength and the free path of electrons. The free 
path depends on the pressure and character of the gas. 
The calculation from first principles of all this process is 
not an easy task. Fortunately, one can measure it quite 
simply by a procedure developed by Townsend.^ We can 
then obtain the effective cross-section ao for a given ex- 
ternal Eq electric field, and use Townsend approximation 
ao a{\£\/Eo) = ao exp{--Eo/\8\) to write the following 
balance equations for electrons and ions 



drUe -I- Vr • je = |Mef"-e| ^0 a{Eo/\£\) (1) 

drUi + Vr • ji = \fieSne\ a{Eo/\£\) (2) 

where ne,i and je.i are particle densities and currents of 
electrons and ions, respectively, and £ is the electric field. 
The fact that the source terms at the right hand side of 
the equations are equal is due to charge conservation in 
an ionization event. 

(m) The electron particle current je is approximated as 
the sum of a drift and diffusion term 



(3) 



where /ig and Z?e are the mobility and diffusion coefhcient 
of the electrons. For anode-directed streamers the ion 
current can be neglected because it is more than two 
orders of magnitude smaller that electronic one, so we 
will take 



0. 



(4) 



{in) The modification of the externally applied electric 
field through the space charges of the particles according 
to the Poisson equation 



Vr • £ = e(n,j - 7ie)/eo- 



(5) 



It is this coupling between space charges and electric field 
which makes the problem nonlinear. 

We want to add a few remarks. In the source term, 
ionization due to the photons created in recombination 
or scattering evens is neglected. This can be justified if 
the cross section of the photoionization process is much 
smaller than that due to electrons. Note that photoion- 
ization can be taken into account, but the dynamical 
equations will become nonlocal. In attaching gases like 
oxygen, a third kind of charged species needs to be taken 
into account, namely negative ions formed by a neutral 



molecule catching a free electron.--- The equations are 
deterministic and stochastic effects are not accounted for 
in the model. 

Finally, the model must be complemented with appro- 
priate boundary and initial conditions. Boundary condi- 
tions will be discussed in detail in the next section. For 
initial conditions, we ignore details of the plasma nucle- 
ation event (e.g. triggering by radiation from an external 
source), and assume that at i = a small well-localized 
ionization seed is present. We also make it clearer below. 

In order to identify the physical scales and the intrinsic 
parameters of the model, it is convenient to reduce the 
equations to dimensionless form. The natural units of the 
model are given by the ionization length i?o = ct^^ , the 
characteristic impact ionization field Eq and the electron 
mobility /ig, determining the velocity vq = figEo and the 
time scale tq = Rq/vq. The values of those quantities for 
nitrogen at normal conditions are 

ag ^ « 2.3 fim, Eq « 200 kV/m, /ig « 380 cmVVs. 

Hence we introduce the dimensionless coordinates^ r = 
R/i?o and t = t/tq, the dimensionless field E = £/Eq, 
the dimensionless electron and ion particle densities a = 
Ue/no and p = Ui/no with uq — eo-Eo/(ei?o), and the 
dimensionless diffusion constant D = Dg/ (RqVq). 
After this rescaling, the model reads 



dta - W-j - a/(|E|) 

dtp ^ a/(|E|) 

p- cr = V • E , 

cr E + D V(T = j . 



(6) 
(7) 
(8) 
(9) 



The function /(|E|) due to Townsend's expression yields 
/(|E|) = |E| a(|E|) = |E|exp(-l/|E|) (10) 

III. NUMERICAL SIMULATIONS 

In this section we will present details of numerical sim- 
ulations of the streamer model discussed previously. In 
confined geometries streamers usually have a nontrivial 
finger like shape. In general two regions can be ob- 
served. The interior of the streamer is an ionized region, 
quasineutral and equipotential. The outer region is filled 
with the nonionized gas. Those two regions are separated 
by a very narrow region in which all the most of the ion- 
ization process is taking place. In this same space there 
is a nonzero charge density and consequently a very large 
electric field gradient. This is one of the reasons why ac- 
curate simulations are rather demanding. These features 
are strongly reminiscent of what occurs in combustion 
frontsii and viscous fingeringii^ 

Simulations based on this model as far we know 
were accomplished by Dhali and Williansi^ and by 
Vitello et. al.^'^ There is also some work by Wang and 
Kunhardt.^^ In Fig.^we can see some simulations of the 
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t = 300 



t = 365 



t = 420 



t = 450 




FIG. 1: Evolution of spontaneous branching of anode directed streamers in a strong homogeneous background field at times 
t = 300, 365, 420 and 450. Model, initial and boundary conditions are discussed in the text. The planar cathode is located at 
z = Q and the planar anode at z = 2000 (shown is < z < 1400). The radial coordinate extends from the origin up to r = 2000 
(shown is < r < 600). The lines denote levels of equal electron density a with increments of 0.2 as indicated by the labels. 



model. A planar cathode is located at z = and a planar 
anode at z = 2000. The stationary potential difference 
between the electrodes A$ = 1000 corresponds to a uni- 
form background field E = —0.5 ez in the z direction. 
For nitrogen under normal conditions, this corresponds 
to an electrode separation of 5 mm and a potential differ- 
ence of 50 kV. The unit of time tq is 3 ps, and the unit of 
field Eq is 200 kV/cm. We use D = 0.1 which is appropri- 
ate for nitrogen, and assume cylindrical symmetry. The 
radial coordinate extends from the origin up to r = 2000 
to avoid lateral boundary effects on the field configura- 
tion. As initial condition, we used an electrically neutral 
Gaussian ionization seed on the cathode 

air, z,t^O)= pir, z.t^O)^ lO"'^ ^-(z-+r-)/im^ 

(11) 

The parameters of our numerical experiment are es- 
sentially the same as in the earlier simulations of Vitello 
et ali^, except that our background electric field is twice 
as high; the earlier work had 25 kV applied over a gap of 
5 mm. This corresponded to a dimensionless background 
field of 0.25, and branching was not observed. Further 
details on this simulation can be found in the referencesi^ 

There have been some more simulations with improved 
numerics and different boundary and initial conditionsiSi 
where this branching phenomena has been observed. 



IV. THE ONE DIMENSIONAL STREAMER 
EQUATIONS IN A COMOVING FRAME: THE 
PLANAR FRONT 

In the previous section some numerical evidence of 
branching have been presented. In this section and the 
ones to follow we will investigate this issue by analyt- 
ical means. Here we will start with the solution for a 
stationary planar front. The idea is to find a uniformly 
translating front and investigate how transversal pertur- 
bation of this solution will develop. 



For planar fronts, we assume that charge varies only 
in the z direction, so using the equations ((BJ-© we can 
write 

dta-dz{aE) - Dd^,a-af{\E\) = , 
dtp - c7f{\E\) = 0, 
dzE - p + a = . (12) 

Next we will change our reference frame to a frame 
moving with velocity v in the z direction (x, y,^ ~ z — vt). 
Then equations (|12|) read 

dta^vd^a + d^{aE) + Ddlcj (jf{\E\), 
dtp = vd^p + af{\E\), 

d^E-p + (7 = 0. (13) 

A front translating uniformly with velocity v in the 
fixed frame is stationary in this comoving frame, dtcr = 
dtp = 0. As a result, the corresponding front profiles are 
solutions of ordinary differential equations. 

We need to set the boundary conditions. The field, 
being completely screened in the ionized region, is ap- 
proximately constant in space and time far ahead of the 
front, so it follows 

E=<f° z^-oo 

\ Eoo Z , Eoo <0 z^+oo ' ^ ' 

where z is the unit vector in z direction. These bound- 
ary conditions imply, that a time independent amount of 
charge is travelling within the front, and no currents flow 
far behind the front in the ionized regime. 

Now, for any nonvanishing far field Eoo, there 
is a continuous family of uniformly translating front 
solutionSfSiii since the front propagates into an unsta- 
ble state»i^ In particular, for Eoo > there is a solution 
for any velocity u > 0, and for E^o < 0, there is a so- 
lution for any v > \Eoo\- These solutions are associated 
with an exponentially decaying electron density profile: 
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an electron profile that asymptotically for large ^ decays 
like cr(^) (X e-'*'^ with A > 0. 

It will "pull" an ionization front along with the same 
speed. (For Eao > 0, the same equation applies for all 
A > f{Eoo)/Eoc, hence for v > 0). For the interested 
reader we refer him to the bibliographypi^ 

Dynamically, the velocity is selected by the initial elec- 
tron profile&iS If initially the electron density strictly 
vanishes beyond a certain point (corresponding to 
A = oo above) 







for C > Co at t = 0, 



(15) 



then this will stay true for all times t > in a coordinate 
system moving with velocity v ~ \Eao\, and an ionization 
front propagating precisely with the electron drift veloc- 
ity |i?oo| develops. In the remainder of the paper, we will 
consider this particular case. 
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FIG. 2: Electron density a, ion density p and electric field E 
for a negative ionization front moving with v = |i5oo| in the 
comoving frame and. The far field is Eoo ~ ~i and D = 0.1. 

In Fig. 121 we have solved H13|) with the boundary con- 
ditions discussed previously (|14|l and (|15|l . We have cho- 
sen the far field Eoc = — 1 and the diffusion coefhcient 
D — 0.1. It was done using a shooting method for solving 
our two point boundary value problem. This technique 
consists in choosing values for all of the dependent vari- 
ables at one boundary. These values must be consistent 
with any boundary conditions for that boundary, but oth- 
erwise are arranged to depend on arbitrary free parame- 
ters whose values we initially randomly guess. We then 
integrate the ODEs by initial value methods, arriving at 
the other boundary. In general, we find discrepancies 
from the desired boundary values there. Now we adjust 
the free parameters at the starting point that zeros the 
discrepancies at the other boundary. The idea is to iter- 
ate this procedure until obtaining the desired accuracy. 
There is a nice pedestrian explanation of solving bound- 
ary value problems by shooting in Numerical Recipes. 



V. SHOCK FRONTS 

In this section we will simplify a bit more our model 
by taking the limit — > 0. For negative fronts, the 
limit Z? ^ is smooth and eliminates the algebraic 
relaxatiouii^ It also reduces the order of the equations. 
We therefore make D = in the streamer equations. 
Then, in the comoving frame, using (|13|l we can write for 
a stationary front 



vd^p 
dfE 



d^{aE) 
<yfi\E\) 

P + (T ^ 



= 0, 



0, 







(16) 
(17) 
(18) 



We can solve this system analytically. If we take the 
(|17|l and subtract it from (|16|l . using H18|l to eliminate 

a — p, we get 



vd^E + aE ^0. 



(19) 



This equation is just a consequence of the charge conser- 
vation. We can see this by writing 9tg-|- V-jtot = 0, with 
the total charge defined as q = p — a. In our model, each 
ionizing collision, produces the same number of negative 
and positive charge, so we end with V-jtot = 0. The total 
current is given by jtot = 9tE-|-fTE and for a planar front 
with constant and time independent field E = EoqZ (|13f) 
in the non-ionized region where ct = 0, the total current 
jtot = jtot{t)z vanishes. In the comoving frame of Eqs. 
(ini and (Iini-lIHll, this means 

The front equations now reduce to two ordinary differ- 
ential equations for cr and E 

d^iiv + E)a] = -<jJ{E) , J{E) = \E\a{E) , 
vd^\n\E\ = cr, (20) 

that can be solved analytically to give 



a[E] 
P[E] 

6-^1 



V 



P[E], 



v + E 
'^-l/(. 



(21) 



dx 



X 



aix)dx, (22) 



v + x dx 



(23) 



This gives us a and p as functions of E, and the space 
dependence E ~ E{(_) implicitly as ^ = ^{E) in the last 
equation. 

We have plotted in Fig.Olthe solutions l(2U-(I2SIl for a 

shock front moving with v = 1. We have chosen ^1=0 
and then -E(Ci) = E^o 



VI. STUDY OF INSTABILITIES: 
CORRUGATION OF THE FRONT 

In this section we will study the stability of the planar 
shock front. It may be unstable with respect to pertur- 
bations having periodicity on the surface of discontinuity 
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FIG. 3: Electron density cr, ion density p and electric field E 
for a negative ionization shock front moving with v = |i5oo 
in the comoving frame. The far field is -Eoo = — 1. 



and then forming "ripples" or "corrugations" on that sur- 
face. In that case, we will be interested in obtaining the 
dispersion relation curve to find which mode will grow 
faster and eventually determine the streamer character- 
istic shape. Here we will derive the perturbed equations 
and the boundary conditions. 

Let the planar shock front which propagates into the 
z direction receive a slight perturbation having an arbi- 
trary dependence on the transversal coordinates x and y. 
Within linear perturbation theory, they can be decom- 
posed into Fourier modes. Therefore we need the growth 
rate s{k) of an arbitrary transversal Fourier mode to pre- 
dict the evolution of an arbitrary perturbation. Because 
of isotropy within the transversal {x, 2/)-plane, we can re- 
strict the analysis to Fourier modes in the x direction, so 
we study linear perturbations oc exp(si -I- ikx). The no- 
tation anticipates the exponential growth of such modes. 
Any perturbation will also lead to a perturbation of the 
position of the ionization shock front. So we will intro- 
duce the new variable C, = £, — t ex.p(ikx + st) and the 
ansatz 



a{xX,t) 
p{x,C,t) 
(t){x,C,t) 



Po(C)+epi(C) 



(24) 



where ctq: Po and 00 are the electron density, ion density 
and electric potential of the planar ionization shock front 
obtained in the previous section. Note, however, that 
these planar solutions arc shifted to the position of the 
perturbed front. Substitution of these expressions into 
(|13ll (with D = 0) gives to leading order in the small 
parameter e 

{v + Eq) d^(Ji = {s + 2(7o - Po- f) 0-1 

-(To Pi + (9c(To - (To/') d^(f)i - sdi^ao, 
V dc_pi = -f <Ji + s pi - (To/' 5f (/>! - sdi^po, 



- fc2) 01 = (Tl - pi + k'^Eo. 

In equations (|25|l we denote / 

diEifm 



(25) 



fiEo), /' 



Eq 



and Eq ~ ^ (?f(/'o(C) as the electric field 

of the uniformly translating front. In the third equation, 
the term k'^(/)i comes as a consequence of the dependence 
of the electric potential with x. 

These equations can be written in matrix form as 





sdc<7o/{v + E) 
sd^po/v 
-Ek^ 




(26) 





1- 2(To - / - Po 


-(To 


dc<To - (To/' 






v + E 


V + E 


V + E 




zl 


s 


-<Jof' 







V 


V 


V 




1 


-1 















1 


/ 










(27) 



Note we have introduce an auxiliary field "01 = '^C'^i 
which coincides with the correction for the electric field 
sign reversed to order e. 

Having obtained the linear order perturbation equa- 
tions, we are now in position to discuss boundary con- 
ditions. First we consider the boundary conditions at 

— 0. There are two types of boundary conditions, 
some arising from the boundedness of densities to the 
left of the shock front at C T 0, and some arising from 
the continuity of fields across the position ^ = of the 
shock front. From 1)16(1 we gather that {v + E) dz<J 
is finite for all z, also for z ] Q and for z = 0, since 
{v + E) dzcr — a {a — p — f) is finite. The same is true 

for (To. In particular, J^^ dz {v + E) dzf^o ^ as Z ^ 0, 
and [v + E) dz<Jo ^ as z t 0. 

Therefore we impose the same conditions for cri, 
namely 



lim 

1^0 



dC {v + E) d^di ^ 



^lim {v + E) d^cTi^O 



(28) 
(29) 



In a second step we are going to make use of the con- 
tinuity conditions. We match the ( < solution to 
the C > solution. As in front of the shock there are 
not sources, one has to solve = for C > and 

V(/) = —Eao z = V z when C, oo. The solution to first 
order in e has the form 



(T = 
p = 



for C > (30) 
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with the undetermined integration constants a and h. 

Now p and V(/) have to be continuous across the shock 
front: V(/) is continuous because the charge density p — a 
is finite everywhere. The continuity of p we get from (|17f) 
and the fact, that a and |E| arc bounded for all z. 

From the continuity of p turns out 

lun(p(a:,C+,i)-p(x,r,t)) = 

^ pi(0) = (31) 

where we have use (|30|) and (|24|l to the right and left 
limits. 

The continuity of the electric field to first order in e 
implies that 



lim 



(ac</)(x,C,t)lc+ -ac0(x,C,t)Ic-) =0 

i(a,0(x,C,<)lc+ -a.0(x,C,i)lc-) -0 (32) 



lim 

c- 

Using expressions (|3U|I and H24|l again, these conditions 
turn out 



(33) 



If we impose the continuity of the potential, we get 
a = (/'o(O) and (/>i(0) = w + b (which is the same condition 
obtained from the continuity of the electric field). 

Finally, from H28|) and H29|) . and taking into account 
that / = (To when C 0, we have 



V'l(0) = S , (7l(0) = 



sff 



(34) 



Collecting all the identities H31() . (|33|l and (|34(l we get 
for the limit of ^ T 




sfiv)/il + s/fiv)) 

z^Q I 

{vk — s)/k 



(35) 



The other boundary conditions, at ( — — oo are the 
total charge equals to zero and the electric field vanishes, 
so they read 



(36) 



where tr^ and (l)^ are constants. 

VII. DISPERSION CURVE 

In the preceding section we have formulated an eigen- 
value problem. Given fc, we want to find s{k) such that 







(^1 


\ 


pi 




f^l 













01 / 









we can find a solution for the transversal perturbation 
equations fulfilling the boundary conditions derived 
previously IpT^ and lp!H|. In general, an analytic treat- 
ment for any value of k is not possible and one has to 
resort to numerical calculations.^*^ However, in the limits 
of small and large wave number the equations simplify 
and we can obtain the asymptotic behaviour of the dis- 
persion relation s{k). 

We will start looking at the small fc-limit. If expres- 
sions (|26|) and (|27|l are evaluated only up to first order 
in fc, then 01 decouples, and we get 




dr pi 



where 



dccj/{v + E) \ 

dcp/v +0(fc2) , 

J 

(37) 





/ s + 2a-f-p 


-(T 


d^a - af \ 






V + E V 


+ E 


v + E 








s 




+ 0{k') 




V 


V 


V 






[ 1 


-1 














(38) 


is the truncated matrix M^^fc 




The fourth decoupled 


equation reads 










9c 01 = 


= Vi 




(39) 



The boundary condition (|35|l turns out 

+ 0(fc2) (40) 



and 









CTO ^ 








01 (0) = 




vk 



(41) 



The expressions (j39|l and (|41|l give a condition on 

■01 



vk 



MO dc. 



(42) 



Consider now the limit s <C f{v). Then Eqs. H37|) 
and H40() up to order s/ f{v) become identical to the per- 
turbed equations obtained from an infinitesimal change 
of Eao- If we compare two uniformly translating fronts 
with infinitesimally different field Eoo at identical posi- 
tions, their linearised difference solves the same equa- 
tions. In this case, "01 is independent of s and k. But 
then H42f) implies 



vk + 0{k^) for fc < Q!(u) 



(43) 
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This result also has an immediate physical interpre- 
tation: 1/fc is the largest length scale involved. It is 
much larger than the thickness of the screening charge 
layer. Therefore the charge layer can be contracted to a 
(5-function contribution along an interface line. Such a 
screening charged interface precisely has the instability 
mode s = vk. 

In the opposite limit, when k becomes large enough, we 
can also find a relation for the dispersion curve. We will 
need to make the assumption that the ion and electron 
densities remain bounded. Taking this into account, we 
can write using (|26|l the equations for t/^i and (jji as 



k'i 

^1 



(44) 



On the short length scale 1/fc, the unperturbed electric 
field for C < can be approximated making an asymp- 
totic expansion of by^" 

E^^v-f{v)C, (45) 
Inserting this expression in 144(1 . we obtain 



0i = fc'(0i-"-/(f)c). 



(46) 



The boundary condition fixes (f>i{0) — {vk ~ s)/k 
and "01(0) — dc^(t)i = s. The unique solution of H4()|) with 
these initial conditions is 



MO 



-feC 



(47) 



for C < 0. Now the mode e^^^ would increase rapidly 
towards decreasing C, create diverging electric fields in 
the ionized region and could not be balanced by any other 
terms in the equations. Therefore it has to be absent. 
The demand that its coefficient (/(w) — 2s)/2fc vanishes, 
fixes the dispersion relation 



for k^ a{v) 



(48) 



Again there is a simple physical interpretation of this 
growth rate. The electric field can be approximated in 
leading order by 



^ ' ^' ' [-2 1; for C > 



(49) 



When the discontinuity propagates with the local field 
V = —E, a perturbation in a field E = —z(v + d^E C) 



will grow with rate dc^E. The averaged slope of the field 
for C > and C < is d^E = f{v)/2, and this is precisely 
the growth rate 1)48(1 determined above. 

We have studied the (in)stability of planar negative 
ionization fronts against linear perturbations and we have 
found 



m = 



\Eoo\ k 

\Eoo\ a{\E^\)/2 



for fc < a{\Eo 
for fc > a(\Eo 



(50) 



So the planar front becomes unstable with a linear 
growth rate s(fc) for small fc to a saturation value 
|i?oo| Ck(|-Eoo|)/2. This gives us a mechanism for branch- 
ing. In the case of a curved front, if the radius of cur- 
vature increases, the planar approximation for the tip 
is sensible and allows a qualitative understanding of the 
branching phenomena. 



VIII. SUMMARY AND OUTLOOK 

In this paper a fully deterministic model for stream- 
ers, without photoionization, which is suitable for nonat- 
taching gases like nitrogen has been presented. We have 
proposed that an anode directed front can branch sponta- 
neously according to this model due to Laplacian interfa- 
cial instability. We have shown some numerical evidence 
of this phenomena. We have studied the stability of a pla- 
nar front and how transversal perturbation would grow. 
This gives us a qualitative picture of the mechanism act- 
ing on a curved front, and we have got the asymptotic 
behaviour of the dispersion curve. 

However, some questions remain to be answered. From 
the dispersion curve any short enough wave length in- 
stability will grow. We do expect that a regularization 
mechanism should come into play. This regularization 
mechanism which selects a particular mode could be the 
electric screening due to curvature. In the present it is 
under investigation.^" Other possibility could be the dif- 
fusion phenomena not considered in the shock front case. 
Diffusion was neglected to prevent mathematical chal- 
lenges, but soon or later one has to face challenges. 

In any case, the physics of low temperature plasmas 
is an area where many fundamental questions are still 
open, where ideas from patter formation, electrodynam- 
ics, quantum mechanics, statistical mechanics and non- 
linear mathematics can be applied, and where the exper- 
imental side has been ahead of the theoretical one. When 
I hear some pessimistic voices for the future of physics, I 
always think there is much room at the bottom... 
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